EIGENVALUES  OF  RANDOM  MATRIC  i'U>  PITTSBURGH  UNIV  PA 
CENTER  FOR  MULTIVARIATE  ANALVSIS  Z  D  BA I  ET  AL  JUN  87 
UNCLASSIFIED  TR-87-16  AFOSR-TR-87-1076  F49628-85-C-8B88  F/G  12/2 


,*IC*OCOPY  RESOLimOfN  TEST  CKWT  l 

Of  STM««Wy2a±-^  ' 


AD-A186  387 


OIK.Eiii  CQEi 


AFOSR.TR.  8  7-  1076 


)N  THE  ASYMPTOTIC  JOINT  DISTRIBUTIONS  OF  THE 
EIGENVALUES  OF  RANDOM  MATRICES  WHICH  ARISE 
UNDER  COMPONENTS  OF  COVARIANCE  MODEL 

Z.  D.  Bal,  P.  R.  Krishnaiah  and  L.  C.  Zhao 


Center  for  Multi variate  Analysis 
University  of  Pittsburgh 


Center  for  Multivariate  Analysis 


University  of  Pittsburgh 


ON  THE  ASYMPTOTIC  JOINT  DISTRIBUTIONS  OF  THE 
EIGENVALUES  OF  RANDOM  MATRICES  WHICH  ARISE 
UNDER  COMPONENTS  OF  COVARIANCE  MODEL 

Z.  D.  Bal,  P.  R.  Krlshnalah  and  L.  C.  Zhao 


Center  for  Multivariate  Analysis 
University  of  Pittsburgh 


June  1987 


Technical  Report  No.  87-16 


\  ^ 


Center  for  Multivariate  Analysis 
Fifth  Floor  Thackeray  Hall 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 


DTIC 

ELECTE 


Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research  (AFSC). 
under  Contract  F49620-85-C-0008.  The  United  States  Government  Is 
authorized  to  reproduce  and  distribute  reprints  for  governmental  purposes 
notwithstanding  any  copyright  notation  hereon. 


t W 


OH  THE  ASYMPTOTIC  OOIHT  HISiniP'Ml'JlIS  OP  THE  EIOTYALUES 
OF  RANDOM  MATRICES  *TIKU  ARI.T  "TIOER  COMPWIEHI3* 

OF  COVARIANCE  MlULL 


Z.  D.  Bal ,  P.  R.  Ktlshnalali  and  L.  C.  Zltao 


ADS TRACT 


In  this  paper,  the  authors  derived  asymptotic  joint  distributions  of 
the  eigenvalues  of  some  random  matrices  which  arise  under  components  of 
covariance  model . 


AMS  1980  subject  classifications.  Primary  62EZ0;  Secondary  62H15,  62H12. 

Key  Wordsyid  Phrases;  Component  of  covariance  model,  elgenstructure  analysl 


i  : 

-  it'i 

limiting  distribution,  random  matrix. 

U  y* 

0  n 

•  1  J *  ■* 

, 

V 

*3 

r*  * 


Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research  (AFSC). 
under  Contract  F49620-85-C-0008.  The  United  States  Government  Is 
authorized  to  reproduce  and  distribute  reprints  for  governmental  purposes 
notwithstanding  any  copyright  notation  hereon. 


.  UNCLASSIFIED 


UCttWIT  UMilMI  A 1 IOM  W*  (Nit  PAGl  (*hi,  Simla  hilmimii} 


REPORT  DOCUMENTATION  PAGE  \ 


HEAD  INSTRUCTIONS 
HfcKQKK  COMPLETING  t-OKM 


1.  NKCtPlCMT'S  CATALOC  HUMIR 


4;  litU  imtt  * 

On  the  asymptotic  joint  distributions  of  the  1 
eigenvalues  of  random  matrices  which  arise  undep  _ 
components  of  covariance  model 


-  June  1987  JS.z/tyjJ 


87-16 


■ 


.  Bal,  P..  R.  Krlshnalah  and  L.  C.  Zhao 


t.  4IMOWWS  OASAMIATION  NMII 

Center  for  Multivariate 


ti 

ysis 


515  Thackeray  Hall 

University  of  Pittsburgh,  Pittsburgh,  PA  15260 


I*.  MA04T  OATS 

June  1987 


It.  NUMB* «  OP  packs 


•  ».  CMTMkUNl  OPPlCC  NAM  I  AN0  AOOMMI 

Air  Force  Offlc'e  of  Scientific  Research 
Department  of  the  Air  Force 

it.  d>»  nr  OM“V9  I  lfi 


14.  MONiromlte  ACCnCV  NAMC  4  aMMIVH  «/(n«U  Kmi  CmMIIMJ  Olht»)  u.  UCUftrT  CLASS.  (»(  |A l,  ifpmfll 

I  \  Unclassified 


Olt'CL  Attlinc  ATlON^OONNCNAOlNC 
SCMKOULC 


'  (•/  All  *OKl| 


Approved  for  public  release;  distribution  unlimited 


I?.  OUTMItytiOM  ttATftMlNT  |*f  Ml#  abiiiacl  MiMed  In  ileal  II,  II  mH—m*  tom  i«NN| 


<t.Y  iOROt  (CanOnue  mi  mvmm  HH  II  NNMNrr  m%4  Him illr  if  Wwi  Ru»*er| 

Component  of  covariance  model,  elgenstructure  analysis,  limiting 
distribution,  random  matrix. 


10  AASI N AC  I  (Cmiilimm  an  t(A  SI  Mc«»twr  M<  *r  Nt«* 

In  this  paper,  the  authors  derived  asymptotic  joint  distributions  of  the 
eigenvalues  of  some  random  matrices  which  arise  under  components  of 
covariance  model. 


ON  THE  ASYMPTOTIC  JOINT  DISTRIBUTIONS  OF  THE  EIGENVALUES 
OF  RANDOM  MATRICES  WHICH  ARISE  UNDER  COMPONENTS 
OF  COVARIANCE  MOOEL 


Z.  D.  Bal,  P.  R.  Krlshnalah  and  L.  C.  Zhao 


i 

ABSTRACT 


In  this  paper,  the  authors  derived  asymptotic  joint  distributions  of 
the  eigenvalues  of  some  random  matrices  which  arise  under  components  of 


AMS  1980  subject  classifications.  Primary  62E20;  Secondary  62H15,  62H12. 

Key  Words  and  Phrases:  Component  of  covariance  model,  elgcnstructure  analysis, 

limiting  distribution,  random  matrix. 


OH  THE  ASYMPTOTIC  JOINT  DISTRIBUTIONS  OF  THE  EIGENVALUES 
OF  RANOOM  MATRICES  WHICH  ARISE  UNDER  COMPONENTS 
OF  COVARIANCE  MODEL 


Z.  D.  Bai,  P.  R.  Krishna la h  and  L.  C.  Zhao 

1.  INTRODUCTION 

Consider  the  following  components  of  covariance  model 

21  +  ilj 

for  j«l ,2,... .m^  and  1*l,2,...,k,  where  X^:  p*l,  Is  the  j-th  observation 
of  the  1-th  group,  y  Is  an  unknown  general  mean  vector,  Is  the  random 
effect  vector  of  the  1-th  group,  and  Is  the  error  vector.  It  Is 

•  i 

assumed  that  ojS  are  lid.,  e^'s  are  lid.,  o^s  and  j^'s  are  Independent, 
and  that 

E(?i)  *  *  2» 

Etajaj)  •  f,  E(c||c||)  ■  Ij  *  0,  (1.2) 

*  E^n5u^2  *  *• 

where  v  Is  a  p*p  non-negative  definite  matrix.  In  contrast  to  the  HANOVA 
model  with  the  non-random  effects,  which  Is  called  the  model  I,  the  model 
(1.1)  Is  also  called  the  model  II.  If  mj  ■  mg  »  ...  ■  m^  ■  m,  the  model 

Is  balanced,  otherwise  unbalanced. 

Recently,  the  components  of  covariance  model  has  received  considerable 
attention  In  literature.  Its  Importance  not  only  consists  In  the 
generality  of  the  model,  but  also  Its  connection  to  linear  structural 
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relationship  models  and  the  signal  processing.  For  a  discussion  on 
applications  In  linear  structural  relations,  we  refer  to  Anderson  (1984) 
whereas  we  refer  to  Zhao,  Krlshnalah  and  Bal  (1986b)  for  applications 
In  signal  processing. 

Many  papers  were  devoted  to  the  balanced  components  of  covariance 

models.  Anderson,  Anderson  and  01  kin  (1986)  studied  the  maximum  likelihood 

estimators  (H.E)  of  p,  v  and  z.  under  the  condition  that  4>  >  0  and  Is  of 

•  1  — 

a  maximum  rank,  and  obtained  the  likelihood  ratio  criterion  (LRC)  for 

testing  the  rank  hypothesis  on  v.  Some  of  these  have  been  obtained  by 

Theobald  (1975),  and  Schott  and  Saw  (1984).  Rao  (1983)  considered  some 

test  procedures  for  testing  relationships  between  two  covariance  matrices 

and  these  procedures  are  useful  In  finding  the  rank  of  the  component  of 

covariance.  Some  results  for  the  case  of  unrestricted  covariance  matrix 

were  announced  by  Anderson  (1984),  where  an  extensive  list  of  references 

Is  given.  Some  authors  put  the  determination  of  the  rank  of  'Filn  the 

framework  of  model  selection;  some  references  of  this  aspect  are  given  In 

Zhao,  Krlshnalah  and  Bal  (1986  a,  1986  b).  In  these  two  papers,  they 

also  proposed  some  procedures  to  determine  the  multiplicities  of  the 

eigenvalues  of  the  matrix  YEj1  or  f. 

For  the  unbalanced  components  of  covariance  model,  Bal,  Krlshnalah 
and  Zhao  (1987)  proposed  some  consistent  procedures  for  the  determination 
of  the  rank  q  of  the  random  effect  covariance  matrix  One  Important  case 
under  consideration  Is  that  Ej  Is  arbitrary  definite  matrix  and 
0  <_q  <  p.  It  Is  assumed  that 

1  I  1*1,2,...,  (1.3] 


and  that 
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k/N  <  y  <  1, 


K 

where  N  =  T  m . ,  M  and  y  are  constants,  M  an  Integer.  We  shall  keep 
1.-1  1 


this  assumption.  In  this  paper. 


Write 


l  BI  X../m.,  X*  -  J  nil/N, 
-1  jSr1J  1  •  i«i 1-1 


(1.5) 


k  m1 


Wi  ^1{lilj"Il^iflj“Il),» 

M2*  < 

Vlg  and  W1  are  called  the  between  groups  and  within  group  sums  of  squares 
and  cross  products  matrices  respectively.  It  Is  easily  seen  that 


(1.6) 


1CT  E(W1>  *  rl*  TuT  E<V  mh*  ^*"7 


We  define 


r2-Il+¥’- 

which  equals  •jj^tfWg)  when  ml  *  **2  *  •••  “  mk* 


(1.7) 


In  the  above  literatures,  elgenstructure  method  occupies  an  Important 
position.  Especially,  the  study  for  the  elgenstructure  of  the  random  matrix 

T^T  V1  the0*r®t1cal  and  practical  Interest.  In  this  paper, 

we  shall  seek  for  the  limiting  joint  distribution  of  Its  eigenvalues. 
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For  the  model  I,  the  joint  distributions  or  asymptotic  joint  distri¬ 
butions  of  the  multivariate  analysis  of  variance  (MANOVA)  matrices  have  been 
found  by  some  authors.  Fisher  (1939),  Hsu  (1939)  and  Roy  (1939)  have 
Independently  derived  the  joint  distribution  of  the  eigenvalues  of  the 
MANOVA  matrix  In  the  central  case.  Hsu  (1941)  derived  the  above  distribu¬ 
tion  In  the  noncentral  case  when  the  sample  size  tends  to  infinity  and  the 
underlying  distribution  Is  multivariate  normal.  Bal,  Krishnaiah  and  Liang 
(1984)  extended  this  result  to  the  case  when  the  underlying  distribution 
is  not  necessarily  multivariate  normal.  In  proving  these  results,  they 
all  assumed  that  the  ratios  of  the  sample  sizes  of  the  groups  to  the  total 
sample  size  tend  to  constants  In  the  limiting  case,  whereas  the  number  of  the 
groups  Is  fixed.  On  the  contrary.  In  order  to  seek  for  the  limiting  distri¬ 
bution  of  the  eigenvalues  of  the  model  II,  we  need  to  assume  that 

the  number  k  of  the  groups  tends  to  infinity. 

The  balanced  and  unbalanced  cases  are  discussed  in  Sections'.  2  and 
section  3  respectively. 
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2.  THE  BALANCED  CASE 


In  this  section,  we  seek  for  the  asymptotic  joint  distribution  of  the 

u  k  i 

eigenvalues  of  the  matrix  -j^^i  ^or  balanced  case»  Here  we  assume 


that  k  •+  »  arid 


m^  "  ®2  ~  • • •  m^  —  m  ^  2 . 


(2.1) 


This  problem  can  be  reduced  to  the  problem  of  seeking  for  the  asymptotic 

joint  distribution  of  the  eigenvalues  of  SgS^1,  where 

n. 


si  *  ^  &  I,!,1 


n 


(2.2) 

I 

S2*T  <2-3) 

and  Y. ,Y„,...,Y  ,  Z1t...,Z„  are  independent  pxl  vectors,  Ylt...,Y_  have 

-n  -1  -n2 

a  common  distribution,  so  do  2lt...,Zn,  satisfying  the  following  conditions: 

EYj  -  EZj  »  0, 

EIl!l  *  zi  >  °»  EZ1Z1  *  e2  ~  °»  (2.4) 

E(lili)2  •=  *•  E(?ih)2  *  “• 

Without  loss  of  generality,  we  can  assume  that  z^  *  Ip  and  z ^  *  diag 
[Xj,...,xp]  with  \l  ^  ...  >.  Xp  >  0,  where  Xj  >_  ...  >_  ^  are  eigenvalues 
of  tgZj1* 

Write 

x1  *  ^  1  ••  •  •  *®h *  h*l,2,...,H, 

where  Cj  >  ^  >  •••  >  ^  0,  ag»0,  a^  *  a^_j+  uh, 

aH  "  p  “  W1  +  +  V 


.  6 

2  2 

Assume  that  n,  varies  as  n  and  11m  n/n.  =  e  ,  0  <  e  <  «  .  Let 

rv~>°  1 

>_  4>pn^  denote  the  eigenvalues  of  SgSj*.  We  are  Interested  In 

the  asymptotic  joint  distribution  of  1*1,2,. ..,p,  where 

c|n)  -  /n(*(")-xi). 

We  need  the  following  lemmas: 

LEMMA  2.1.  Let  Xn,  n=0,l,2,...,  be  a  sequence  of  random  p-vectors  with 

X„  In  distribution.  Then  there  exists  a  probability  space  (n,F,P) 

on  which  we  can  define  a  sequence  of  random  vectors  X.  n=0,l,2,...»  such 

**n 

that 

1.  Xn  and  Xn  are  identically  distributed. 

2.  +  >L  polntwise. 

n  «u 

The  above  lemma  is  given  in  Skorokhod  (1956). 

LEMMA  2.2.  Let  9n(x)  be  a  sequence  of  K-degree  polynomials  with 

roots  x<n),...,x<">  for  each  n,  and  let  g(x)  be  a  k-degree  polynomial  with 

roots  Xj,...,X|(,k  <_  K.  If  9n(x)  -*■  g(x)  as  n  ®,  then  after  suitable 

rearrangement  of  xjn^,... ,x£n\  we  have  xjn^xj»  j=l»...»k  and 

|x<">|  ->  -,  j*k+l,...,K. 

For  a  proof,  we  refer  to  Bal  (1984). 

-2  -2  -1 
Put  n^  *  e  nBn  with  en  1.  The  eigenvalues  of  SgSj  are  the  roots 

of  the  determinant  equation 


det  (Sg-^Sj)  *  0. 


(2.5) 


Write 


u„  -  <»<J>>  -  ^(srri)  .-jr  j/Wl-Ip). 


Then  (2.5)  becomes 


det  (£  *V»  '  O. 


(2.6) 


where 


0=  /(*-C1)Iy1 

By  the  central  limit  theorem  and  Lemma  2.1,  we  can  assume  that 


Un  ’  U  4  a-s‘ 

V_*  V  4  (v.i).  a.s. 


(2.7) 


n  ’  '/ij 
satisfying  the  following  conditions: 

(1)  {u^,  l^i  <  j  <  P>  and  {v^,  1  <.  i  <_  j  p}  are  two  p(p+l)/2 

dim. normal  vectors  independent  each  other. 

(2)  u1d  -  uJt,  vfj  =  Vj,  for  all  1,  j. 

(3)  EU...  =  Evi j  =  0  for  all  1,  j. 


U 


and 


(4)  covtv^.Vjj)  =  EZ^2  Z.j  -  \^\y  if  1  <  1  IP* 
cov(uii  »ujj)  =  eyh2yi/  -  if  1  li  i  J  IP* 


C0V(V1112.VJ1JZ)  = 


cov(u11i2»uj1j2)  "  EY111Y112Ylj1Ylj2» 


if  1  1  1  i2  -  p*  1  1  1  ^2  -  p  anc*  at  'east  one  of  the  strict 

inequalities  i^  <  i2  and  <  j2  holds.  Here  and  are  the  i-th 
components  of  the  vectors  Yj  and  Zj,  respectively. 

Split  the  matrices  Un,  Vn,  U  and  V  into  blocks  as  follows 

Un  *  <C>-  »„-(»&’*  u-(Uhg),  »•(»„)• 

h,g  -  1  *2 ». . . »H, 

where  ujg\  V^,  Uhg,  Vhg  have  order  uh*ug,  h,g  =  1,...,H. 

Take  the  variable  transformation  <t>  =  ^  +  s//rT.  Multiply  by  n*^ 
the  first  pj  rows  and  the  first  p^  columns  of  the  determinant  on  the  left 
hand  side  of  (2.6),  (2.6)  can  be  rewritten  as 

fn(c)  -  0,  (2 

where  fn(s)  is  a  polynomial  in  x,  with  degree  p,  and 


f  (c)  ■  f(c)  ■ 

fH«o 


(^2"^)^2 


After  the  variable  transformation,  the  p  roots  of  (2.6)  (i.e., 
can  be  written  as  s  »^rr( n^-Cn ) -  Since  £im  4>jn^  =  *.• » 

•  *  i  *  n-*«o 


(2 


(2.8)) 

1=1,2,... 


we  have  £im 


i=a,  +  1 


P 
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By  Lemma  2.2, 


-£1m  cL ^  i 
n-*» 


(2.10) 


where  Ci are  the  roots  of 
l  ax 


det  )  *  0. 


(2.11) 


’P 

are  the 


Similarly,  we  can  prove  that  the  joint  distribution  of  e{n^,..,Cpn^  tends 

to  that  of  as  rv*«,  where  ca  +1  ^  ?a  +2  c 

^  h-1  h-1  “f 

roots  of 

det  (Vhh-6ChUhh-cIUh)  =  0,  h=l,2,...,H  (2.12) 

Thus,  we  have  proved  the  following 

THEOREM  2.1.  Let  S^Sg  be  defined  by  (2.2)  and  (2.3)  respectively  and 
let  <J>|n^,...,<frpn)  be  the  eigenvalues  of  SgSj"1.  Write 
=  *ffi($jn^-Xj),  1=1,2,... ,p.  Suppose  that  Yi»***»Yni»zi»***»zn  are 

independent  and  Y.,...,YM  have  a  common  distribution,  so  do  Z,,...ZM, 

•!  ~n 

where  EYj  =  EZj  =  0,  EY^J  =  Ip,  EZjZj  =  diag  [Xlt...,Xp],  E(yJYj)2  <  »,  E(ZJZ1)2<»I 

and  that  -Elm  —  =  e2  >  0.  Then  the  joint  distribution  of  ?!n^,... ,c*n^ 
n->-»  nl  1  p 

tends  to  that  of  c, ,...,tn,  where  >  ...  >  ca  are  the  roots  of  the 

1  p  ah-l+ 1  “  ~  ah 

determinant  equation 

det  (QLL-cl  . 

V 


det  ^hh'51-  )  ’  °> 


(2.13) 


and 


satisfying 


Qhh  *  ^ij^*  1  »j-a^_j+l ,. . . ,a^,  h*l ,2 ,. . . ,H s 
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0)  qij  =  qj1  for  each  (i,j). 

(2)  {q^j.  ah-i+1  i  j  -ah*  h=1»2 . H}  has  a  joint  normal 

distribution. 

(3)  Eq^.  *  0  for  each  (i,j), 

Cov(qii.qjj)  -  -  5hS(1+e2)- 

if  ahl  +  1  <  1  <  ah,  agl  +  1  <  j  <  ag,  1  1  h  1  9  i  H» 


and 


Cdv(q1112>qj1j2)  ■EYli1Y112  Ylj1Ylj2  +  9  5hcgEZli1Zli2Zlj1Zlj2* 

if  ah_j+l  —  M  —  ^2  -  ah*  ag-l+*  —  ^1  —  ^2  -  ag*  *  —  ^  —  9  —  H  and 
at  least  one  of  the  inequalities  ij  <  i2,  jj  <  j2  holds. 

Note  that  (2.13)  is  equivalent  to  (2.12)  by  the  independence  between 
U  and  V,  and  Theorem  2.1  holds  evidently. 

THEOREM  2.2.  Besides  the  assumptions  of  Theorem  2.1,  we  assume  that 

EV11Y1J  ■  *  !• 

EZliZlj  '  EZ11EZ1J  *  Vj*  1  -  1  ‘  J  -  p> 

and 

EYli1Yli2Ylj1Ylj2  =  EZl1iZl12ZlJlZl-i2  *  ° 

for  all  ij,  12,  jj,  j2  satisfying  that  ah-1  +  1  <.  1j  <.  12  i  ah» 

agl  +  l^j1^j2^ah,  l<h^g^Hand  that  at  least  one  of 

ij  <  12,  jj  <  j2  holds.  Then  the  limiting  vectors  (cj,...,i;al),..., 

(?»  .1 ,...»?_)  are  Independent  each  other, 

aH-l+1  p 


THEOREM  2.3.  Under  the  assumptions  of  Theorem  2.2  and  that 
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EY^  *  3,  EZ^  -  3xZ,  1-1,... ,p,  the  joint  distribution  of 

n|n^  «  ^(♦|n^-5h)/(4(l+eZ)ch),  1  *ah„i+l ••  •  •  *h*l *2 , . . . ,H  tends  to  a 

limiting  distribution  with  density 


VI 

II  ii  »»»,n_  ) i 

i-l  ah-l+1*  ah 


where  D(x^,...,x  )  Is  given  by 


D(x1....,xu)  -  2-“/2(i;i  '^)'1{y<Mv«i)]  «p[-7 


•  J  X1  ix2  i  ...  > 


(7efer  to  Bai,  Krlshnalah  and  Liang,  (1984)). 

To  prove  this  theorem.  It  Is  enough  to  notice  that  ^(l+e2)]"1^2  Q, 


hh 


Is  a  (real)  central  Gaussian  matrix  of  order  wh*uh  and  that  D(xjt...,x  ) 

is  the  joint  density  of  the  eigenvalues  vj  ^  ...  of  a  central  Gaussian 

matrix  of  order  y*u.  For  the  definition  of  Gaussian  matrix,  refer  to 
Zhao,  Krlshnalah  and  Ba1(19S7). 


i 


i 


! 

I 
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3.  THE  UNBALANCED  CASE 

In  this  section,  we  study  the  asymptotic  joint  distribution  of  the 
eigenvalues  of  defined  by  (1.6).  For  simplicity,  we  assume  that 

®1  ejj  In  (1.1)  are  distributed  as  p*l  normal  vectors  for  each 
(1,j),  and  that 

«l's  are  lid.,  ejjS  are  lid.,  c^'s  and  e|jS  are 
Independent,  E(aj)  -  E(en)  -  0,  and  (3:1) 

^(®  1- 1 )  ■f  »  Ekli*ll)  *  Ei  >  o. 

Besides,  we  assume  that 

IPf  *  ®2  ♦  «(~).  6  *  °»  •  constant,  (3.2) 

*  ,1"?  *  "*  +  1  i  »i  i  1.  1*1 . k.  (3.3) 

1-1  Ai  1 

m  and  M  are  constants. 

let  «{k*  >, ...  >  >.  ...  >_  *p  and  xjk^  i  ...  i  x£k)  denote 

the  eigenvalues  of  ^  «*d  yj1  Wectlvely,  where 

h  9  ri  * 1 f  "  h  * 

Without  loss  of  generality,  we  can  assume  that 

Ej  -  Ip  and  t  •  dlag  . *p]. 


xOO  .  *00 

1  h  • 


vv 


If  h»l,2,...,H, 


Write 
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«  •  •  (l\  y  * 

where  *1  >  *2  *  •••  >  '  "  *  +  7  *h*  h«l,2,...  ,H. 


The  eigenvalues  of  ire  the  r00ts  of  the  det«m1nMt 


equation 


det  (t£i*2  -  TCfV  ’  °- 


(3-4) 


Put  N-k  ■  e”2ke’2,  then  -  1  +  o(~)  as  k-*». 
K  *  ✓F 


Write 


uk  ■  '  V* 

vk  -  (v<5>).  ^(^2  -  I2)  *  -  <H»9  . ^k): 


Then  (3.4)  becomes 


d«t(—  »t  -  ^.V,,  -  Dh)  •  0. 
/f  k  Jt  k  * 


0.5) 


where 


“1 


\  -  U  *  (ujj),  a.s. 


0.«) 


where  u^  and  v^  are  all  nora«1  variables  for  all  1,j,  and,  U  and  V  are 


Independent.  Besides,  If  we  split  the  Matrices  Uk,  *k,  U  and  V  Into  the 
fol lowing  blocks 

«k-<«J,k)>.  »k*<Yi.k>)’  U'(V-  v-(V* 

*  1,2,. ..,H, 


with  uj^,  vg>.  U^,  Vhg  being  of  order  uh*vg,  then  Uhh,  Vhh,  h-1,2 . H, 


are  Independent  and,  ^  and  (2m2*jJ  +  4(e’2+l)Jh  +  2)-1^2  Vhh, 

h»l,2,...,H,  are  all  central  Gaussian  matrices. 

Let  ■  «/F(ijk^  -  x<k>).  1-1,2,... ,p.  Using  the  argument  used  In 

the  section  6,  we  can  prove  that  the  joint  distribution  of 

tends  to  that  of  c.,...,t  as  k-*-,  where  c,  >  c,  . ,  >  ...  >  e, 

1  P  *h-i  1  —  ah-l  *  ~  ~  ah 


are  the  roots  of 


det^Vhh"eehUhh“5luh^  *  °»  h-1,2, ...,H, 


(3.7) 


where 


■  1  ♦  0+®  )♦>,.  h-1,2 . H. 

Now  [2m2*2  +  4(e“2+l)*h  +  2  +  2eh]*l^2(VW|-  05hUhh),  h-1,2, ...,H,  are 
Independent  central  Gaussian  matrices.  So  we  have  the  following 


THEOREM 


3.1.  Suppose  that  and  In  (1.1)  are  both  pxl  normal 


vectors  for  each  (1,j),  and  (3.1)-(3.3)  are  satisfied.  Put 


•♦2^2«h. 


1  *Rk_  i  ^1 ,  •  •  •  »a|. ,  h-1,2,...  ,H, 


where  ■  1  ♦  (l+e"2)*^.  Then,  the  Joint  distribution  of  r»j  ,. 
tends  to  a  limiting  distribution  with  density 


' 00 
.  »Hp 


II 

H  0(n_  )i 

h-1  *h-l+1  *h 

where  D(x1,...,xv)  has  given  In  Theorem  2.3. 
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